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Abstract 


The density-matrix renormalization group (DMRG) method, which can deal 
with a large active space composed of tens of orbitals, is nowadays widely used as 
an efficient addition to traditional complete active space (CAS)-based approaches. 
In this paper, we present the DMRG algorithm with a multi-level (ML) control 
of the active space based on chemical intuition-based hierarchical orbital ordering, 
which is called as ML-DMRG with its self-consistent field variant ML-DMRG-SGF. 
Ground and excited state calculations of H 2 O, N 2 , indole, and Cr 2 with comparisons 
to DMRG references using fixed number of kept states (M) illustrate that ML- 
type DMRG calculations can obtain noticeable efficiency gains. It is also shown 
that the orbital re-ordering based on hierarchical multiple active subspaces may be 
beneficial for reducing computational time for not only ML-DMRG calculations but 
also DMRG ones with fixed M values. 
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1 Introduction 


One of the main challenges in the current community of theoretical chemistry is the de¬ 
scription of the systems with complicated electronic structures, e.g. the transition metal 
complexes, open shell and excited electronic states, and the bond-breaking/formation 
00 - In such cases, several leading components exhibit in the wave function due to 
quasi-degeneracy and accordingly the accuracy of the Hartree-Fock (HF) approach and 
those correlation models which rely on the HF assumption are deteriorated. For dealing 
with these quasi-degeneracy problems, multi-reference (MR) approaches [^|^, which use 
a number of important reference states instead of only one in the single-reference (SR) 
approaches as the basis to build further possible excitation conhgurations, are usually 
required. In MR approaches, the orbitals are usually taken from a multi-conhgurational 
(MC) self-consistent held (MCSCF) [5]-[^ calculations. In a most popular MC approach, 
the complete active space (CAS) self-consistent held (CASSCF) method, all possi¬ 
ble conhgurations in the active space are constructed, thus the quasi-degeneracy can be 
taken into account. Thanks to the latest methodological developments 


10-16 , great 


progresses have been made for some of the challenging MR issues such as exploring the 
multi-radical character of one- and two-dimensional graphene nanoribbons (GNRs) 
and analyzing the magnetic interactions in molecules and highly correlated materials . 
However, the size of the active space used in CASSCF is usually no larger than (16e,16o) 
due to the exponential growing of the computational cost with the number of active or¬ 
bitals and electrons [^|^, and this greatly hampers the routine applications of CAS-based 
MR approaches to the large molecular systems. 


In recent years, the ab initio density-matrix renormalization group (DMRG) method 


19 93 emerges as a promising MC quantum chemical approach which can deal with active 


spaces larger than those by traditional CAS-based methods. In the DMRG method which 
originated from the community of condensed matter physics by the pioneering work of 
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White the M eigenvectors with largest eigenvalues of the reduced density-matrix 
(RDM) of sub-systems constitute a truncated basis set for the purpose of reducing the 
freedom of the Hilbert space. Because its computational cost scales only polynomially 
with the number of active orbitals and the number of retained renormalized basis states 
M [^, DMRG is nowadays regarded as an efficient alternative to the full conhguration 
interaction (FCI) or the CAS configuration interaction (CASCI) method which has an 
exponential scaling. With the latest improvements of the DMRG method (e.g. spin- 


adaptation matrix-product state (MBS) [%-98 or tensor network state 

(TNS) [^1^ representation of the wavefunction, and the incorporation of dynamical 
correlations 52,59,^ 77,78,93 as well as the consideration of relativistic effect 33 


86,^), DMRG may be applied to molecular systems with exceedingly large-size CAS 


which are completely unaffordable by traditional approaches. For example, Hachmann 


et ah 39 firstly illustrated the polyradical nature in higher acenes through (50e, 50o) 


DMRG calculations, and later Mizukami, et ah 71 further extended the radical nature 
study to GNRs by (84e, 84o) DMRG calculations and demonstrated that the mesoscopic 
size effect that comes into play in quantum many-body interactions of a large number of 
vr electrons is responsible for the polyradical nature in GNRs. Very recently Kurashige et 


ah 73 presented the hrst ab initio evidence of possible non-adiabatic chemical pathways 
of Mn 4 Ca 05 cluster in photosystem II by virtue of (44e, 35o) DMRG calculations, and 


Sharma et al. 74 critically examined the validity of the consensus phenomenological 
models of iron-sulfur clusters with the help of (30e, 32o) and (54e, 36o) DMRG benchmark 
calculations and found the widely used Heisenberg double exchange model underestimates 
the number of states by one to two orders of magnitude. Nevertheless, it is well-known that 
in many cases the optimal active orbitals cannot be accurately pre-identihed in a simple 
way and accordingly the optimization of active orbital is usually required for MG and MR 
calculations. Therefore, it is of crucial importance to perform DMRG calculations with 
optimized orbitals, which may be self-optimized by energy minimization, along the same 
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lines as CASSCF 44,4^. Alternatively, one could also carefully choose various molecular 


orbitals from other low-level electron-correlation calculations as DMRG basis 55 69,76 


An important feature of DMRG is that it allows to capture all types of electron cor¬ 
relation effects (dynamical and non-dynamical) in a given active space in a consistent 
way. In many cases, choosing a hxed number M of preserved states in DMRG-based 
calculations is feasible but not economic. Legeza et ah suggested a dynamical 

block state selection (BBSS) protocol based on a fixed truncation error ptr of the sub¬ 
system’s reduced density matrix instead of using a hxed number M of preserved states 
in DMRG sweeps. This approach leads to substantial gains in computational efficiency, 
because such a protocol can automatically adapt the computational demand of a DMRG 
step at a specihed orbital to its special electron-correlation environment. Since not every 
renormalized state contributes the same amount to the hnal energy at every point in 
the sweep, a hxed constant M as used in standard calculations can be a bit inefficient 
when dealing with very weakly occupied/unoccupied orbitals. Similarly, the truncation 
by density-matrix threshold similarly is inefficient because a hxed density-matrix weight 
does not mean the same thing energetically at the ends of the sweep as at the middle. 
It is noticeable that the idea of “multi-level (ML) control of the active space” is pop¬ 


ular in the MG approaches 99-106 , such as restricted active space self-consistent held 


(RASSGF) [MTOO and generalized active space (GAS)-GI/SGF [104 -106 . In such kind 
of methods, the whole active space is divided into several subspaces according to their 
distinct electron-correlation environments, and each subspace is controlled at a diherent 
computational level. For example, GASSGF calculations on the Gd 2 , and oxoMn(salen) 


complex by Ma et ah 106 showed that GAS wave functions achieve the same accuracy 


as GAS wave functions on systems that would be prohibitive at the GAS level. However, 
whether the ML concept can be ehectively applied to the DMRG method is still unknown. 


In this paper, we intend to investigate the possibility and the performance of DMRG 
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calculations using hierarchical orbital reordering of multiple active levels based on chemical 
intuitions, as commonly used in RAS calculations. The structure of this paper is as follows. 
In Sec. II, we present a short description of the DMRG and ML-DMRG algorithms as 
well as how to combine ML-DMRG with an orbital-optimization scheme. Illustrative 
numerical examples of the ground and excited states of H 2 O, N 2 , Gr 2 and indole with 
comparisons to standard DMRG references are presented in Sec. IV. Finally, we draw our 
conclusions in Sec.V. 


2 Methodology 


In the MRS ansatz, the DMRG wave function can be expressed as 97 




( 1 ) 




where matrix is for site i with the local state label [cTj], the dimension of 

is M X M at most, and L is the total number of sites. The determination of an MPS 
wave function is a non-linear problem. However, it can be reduced to a linear problem 
by the sweep algorithm, by which the variational parameters in the MPS wave function 
can be variationally determined by just iteratively optimizing one or two of the matrices 


A* while keeping all other matrices hxed 94 . The MPS wave function can provide both 


qualitative and quantitative descriptions of the electronic structure, achieved by using 
small and large M. The polynomial scaling -|- M^L^) of the computational effort 

(in sweeps) with the size of system is owned by DMRG algorithm . Usually the M is 
hxed by a specihc number (or dynamically adjusted in DBSS [^). When applying the 
DMRG algorithm into quantum chemistry, one need to use it to deal with the second 


quantized Hamiltonian 51,62 




( 2 ) 




ijkl 
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where a~^ and a represent creation and annihilation operators, i, j, k, I denote orbitals (e.g. 
HF canonical molecular orbitals) which are more appreciated than “sites” by chemists, h 
and V are the 1- and 2-electron integrals, respectively. 


In our ML-DMRG framework, we further regroup orbitals into different orbital spaces, 
and construct the DMRG wave function for the active space with a multiple M treatment 
as illustrated in Fig.[^ Usually the subspace, which is assigned in the middle of subspaces, 
should contain the most important orbitals for the purpose of maximizing the efficiency 
of the sweep algorithm in DMRG (2^ . 



(T| <7^ ''' fT, • ■ • <T. CT, • • • cr,_| <j j 




Orbital reordering and assortment 

Level-1 orbital group Level-2 orbital group Level-3 orbital group 



fTj. (7^. • • ♦ fr. (7,. • • • C7., (J.. • • ♦ cr,_,. (T^. 


Left -> right: Afi -'Ri ^2 ^2 ^ ^2 ^2 ^ 

Right <- left: M^xM^ M^x x M.^ M3 xM3 


Figure 1: Illustration of ML-DMRG treatment. The dimension of A matrix of the bound¬ 
ary sites(i.e. orbitals) in different sweeping directions are shown at the bottom. 

When combining ML-DMRG with orbital-optimization, a two-step procedure is im¬ 
plemented. At each iteration, an ML-DMRG calculation is performed hrst and a reference 
state is obtained. Then it is followed by orbital-optimization process, according to the su- 
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per configuration interaction (super-CI) scheme 107 . With the new set of orbitals, a new 
ML-DMRG iteration is performed, until there are no orbital rotations, which indicates 
that the ML-DMRG-SCF convergence has been reached. The flowchart of ML-DMRG- 
SGF is also illustrated in Fig. and two types of orbital-rotation treatments are used: 
one is the original treatment, in which all the necessary orbital rotations are allowed; the 
other is the simplihed treatment, in which the orbital rotations between different active 
subspaces are treated as redundant. 



Virtual 


Low-level subspace 


High-level subspace 


Low-level subspace 


frozen/closed 



Figure 2: Treatment of orbital rotations in ML-DMRG-SGF method: 1) Original treat¬ 
ment (left) shows that all the necessary orbital rotations are allowed; 2) Simplihed treat¬ 
ment (right) shows that orbital rotations between different subspaces are treated as re¬ 
dundance. 


3 Computational Details 


We use a locally modihed version of the Block program 63 to implement the ML 


DMRG calculation. When running the ML-DMRG-SGF calculations, we use a develop- 










































ment version of the MoLCAS package 108 , which contains a Molcas-DMRG interface. 
This Molcas-DMRG interface snpports Block and Maqlfis QG-DMRG packages. 
The single point calculations by DMRG-GASGI/SGF, ML-DMRG(-SGF), and geometry 
optimizations can be implemented using this interface. The GASSGF, RASSGF calcula¬ 
tions of PEG of N 2 are implemented by Molpro.2010 [109 , and all the other calculations 
are implemented by the development version of Molcas package. 


4 Results and Discussions 

4.1 Accuracy of ML-DMRG 

Firstly, we investigate the numerical accuracy of the ML-DMRG on the water molecule 
with a ANO-RGG-VDZP basis. In this example, the core orbitals and the valence orbitals 
are chosen as the high-level orbitals (lOe, 7o) and the remaining orbitals are chosen as 
the low-level orbitals (17o). In the ML-DMRG calculations, the number of preserved 
state (M) is hxed to 3000 when sweeping among high-level orbitals and the M value is 
range from 1 to 300 when sweeping among low-level orbitals. For comparisons, we also 
run the GASGI, GASSGF, GASPT2, MRGI calculations as well as coupled-cluster (GG) 
calculations. All the results are illustrated in Fig. It can be seen that the accuracy of 
ML-DMRG is gradually improved with the increasing M value in the low-Ievel subspace. 
If only 1 preserved state is used in low-level subspace, ML-DMRG can reproduce the 
result of GASGI in only the high-level subspace. It means that only the non-dynamical 
correlation in the high-level subspace is considered, while there is no dynamical correlation 
of low-level subspace being involved. When increasing the number of preserved states that 
used in low-level subspace, one would hnd that ground-state energy decreases gradually, 
and hnally the energy value (-76.3267187 Hartree, M=300 in low-level subspace) is even 
lower than that of GGSD(T) (-76.3265256 Hartree), and it is quite close to that of standard 
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DMRG calculation using M=3000 (difference is about 0.0005 Hartree). It means that the 


electronic correlations within the low-level subspace and those between the high-level and 


low-level subspaces can be effectively incorporated upon increasing M used in low-level 


subspace. 



-76.328 


50 100 150 200 250 

M in low-level subspace 


300 



Figure 3: Results of accuracy analysis along with increased number of preserved states M 


in the low-level subspace of the H 2 O molecule. Here 4ai 2bi 0 a 2 I 62 orbitals are chosen as 


the high-level orbitals with M=3000, 7ai 56i 2 a 2 3^2 orbitals are chosen as the low-level 


orbitals with different M (shown in hgure). The ANO-RCC-VDZP basis is used here. 


Secondly, we extend the active space to all the orbitals of N 2 molecule with further 
partition of orbital levels, in order to reveal the capacity and behavior of ML-DMRG 
when dealing with large active spaces (around 60 orbitals). The N 2 molecule in both 
equilibrium geometry (l.ORg, Re=l-0976A) and near dissociation geometry (S.ORe) with 
cc-pVTZ basis are used. The NOs calculated from GASPT2 are used as the basis when 
implementing DMRG or ML-DMRG calculations. In these calculations, the inner or¬ 
bitals {Is), valence orbitals {2s2p), third-shell orbitals {3s3p3d), and forth-shell orbitals 
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{4s4p4d4f) are the orbitals of the four different subspaces. When implementing ML- 
DMRG calculations, the numbers of preserved states in different subspaces are kept by 
the proportion of 0.25:1:0.5:0.25. It is clearly shown that, for moderate energy accuracy 
criterion (no less than 10“^ Hartree) ML-DMRG calculations can have noticeable effi¬ 
ciency gain over standard DMRG calculations with hxed M values (Fig. |^. One may 
also notice that the BBSS strategy with ^^,.=1.0x10“^ has significant computational speed 
advantages for low energy accuracy (around 10“^ Hartree) calculations, which save more 
than 50% computational time. However, more accurate DBSS-DMRG calculations with 
Pir=l-0xl0“^ can be easily out of the capability of current computational resources and 
hardly give converged results because it usually requires too many states to be kept (more 
than 10000, see Fig. Sla ing), Under such circumstances, using fixed M value or the 
ML strategy will be a computational advantage for DMRG calculations. Therefore, one 
can get a useful hint: for achieving computational efficiency gains, the BBSS strategy 
provides a good choice for low-accuracy DMRG calculations, while the ML strategy can 
provide an alternative solution for medium-accuracy DMRG calculations. 


11 







Figure 4: Energy convergence of full orbital DMRG calculation of N 2 molecule at 1.0-Re 
(1.0976 A, left part) and 3.0-Re (3.2928 A, right part) intra-molecule distance. The ML- 
DMRG calculations were implemented with M= 750, 3000, 1500, 750 in four different 
subspaces. Every point of all the four methods represent the difference between the energy 
when sweeping on the middle of active orbitals and the final extrapolated energy. The 
latter parts of these lines are flattening out, showing that the calculations are converged at 
this level. For comparisons, the DMRG calculations were implemented with M=3000 in 
the whole active space as well as two DBSS-DMRG calculations with different truncation 
cutoffs of the reduced density matrix ptr (A: 1.0x10“^; B: 1.0x10“®) are also shown. 

All the calculations were implemented using 16-core Inspur clusters with Intel Xeon E5- 
2670 GPU. 
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4.2 Accuracy of ML-DMRG-SCF 


We also build ML-DMRG-SCF and apply it to the calculation of the N2 molecule. 
Actually, (6e, 60) should already give good descriptions of the dissociation process of N2, 
if only a qualitative picture is required. In ML-DMRG-SCF, the core (is) and valence 
orbitals (2s, 2p) are chosen as the high-level orbitals with (M=2000), while the 3s, 3p, 
3d orbitals are chosen as the low-level orbitals with (M=200). Here we implemented 
two types of orbital optimization strategies as shown in Fig. 1) Orbital rotations 
between different levels are allowed, i.e. original treatment or RASSCF-type treatment, 
and 2) these are not allowed, i.e. simplified treatment. The cc-pVQZ basis is used. 
When the RASSCF-type rotations are turned on, we found that the ML-DMRG-SCF 
may diverge after several SCF iterations. Once the inter-subspace rotations are turned 
off, the convergence behavior became much improved (Fig. |^. In the latter case, the 
results of ML-DMRG-SCF have stable deviations (10 ~ 14 mH) when comparing with 
those by DMRG-SCF (Fig. [^. When the optimized orbitals from ML-DMRG-SCF are 
used as the basis to run further higher-level DMRG calculations, the deviations between 
results of DMRG and those of DMRG-SCF are less than 1 mH. It means the optimized 
orbitals from ML-DMRG-SCF calculation owns a good quality comparable to that by 
DMRG-SCF calculation. It also implies that the ML-DMRG-SCF could be used as the 
alternative to DMRG-CASSCF for yielding optimized orbitals, or be used as a pre-DMRG- 
SCF calculation in order to run a DMRG-SCF calculation more efficiently. 
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Table 1: Calculated energies (a.u.) of N2 system in different geometry configu¬ 
rations by DMRG-SCF and related energy differences (a.u.) between RASSCF“, 
ML-DMRG-SCF^’C DMRG/ML-DMRG-SCF'^ results and them. 


Method 

O.SRe 

l.ORe 

1.3Re 

2.0i?e 

3.0Re 

4.0Re 

RASSCF“ 

0.00824 

0.00866 

0.00239 

0.00823 

0.00849 

0.00856 

ML-DMRG-SCF'’ 

0.0068 

0.0063 

0.0090 

0.0056 

0.0049 

0.0051 

ML-DMRG-SCF= 

0.0112 

0.0132 

0.0138 

0.0106 

0.0101 

0.0101 

DMRG/ML-DMRG-SCF'' 

0.00056 

0.00097 

0.00054 

0.00008 

0.00004 

0.00004 

DMRG-SCF 

-109.13671 

-109.36550 

-109.22427 

-109.01743 

-109.00660 

-109.00618 


^ Orbital rotations between different active subspaces are allowed. Notice that extra attempts (e.g. start from 
pre-optimized orbital) may be used, since it may diverge during the ML-DMRG-SCF iterations in this case. 

^ No orbital rotations between different active subspaces. 

^ RASSCF(14e,28o) with CAS(6e,6o) as RAS2, two holes are allowed in RASl and two electrons are allowed in 
RAS3. 

DMRG(14e,28o) calculations (M=2000) with optimized orbitals from ML-DMRG-SCF^. 
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Figure 5: The convergence behavior of ML-DMRG-SGF (without orbital rotations be¬ 


tween different active subspaces) and DMRG-SGF. 


14 



































4.3 State-averaged ML-DMRG-SCF calculation of absorption 


spectra 


Hereby, we implement the state-averaged ML-DMRG-SCF method in order to evaluate 
the electronic excitations in complicated molecular systems. As an example we choose 
to study the indole molecule (see Fig. |^. The indole chromophore is responsible for the 
low-energy absorption and emission spectra in order to understand the tryptophan spec¬ 


troscopic properties 111, 112 , which are believed to be useful probes of local environment 


and dynamics in proteins 113 . The studies on the electronic spectrum of indole have 
been almost exclusively focused on the position and nature of the two low-lying valence 


excited singlet states {2^A and 3^A'), which were early labeled and states 114 


In a previous study, Serrano-Andres and Roos pointed out that a valence basis set and 
a valence active space do not guarantee that the obtained states have a pure valence 
character, and as such, the recommended procedure is therefore to always include diffuse 


functions in the basis set and Rydberg-type orbitals in the active space 115 



Figure 6: Indole molecule, the geometry was optimized by BP86-D3/augTZ 


State-averaged ML-DMRG-SCF and state-averaged DMRG-SCF are used to calculate 
the excitation energies of indole molecule. When using ML-DMRG-SCF, the choice of 
active subspaces is as following: 1) valence orbitals (9 tt) plus some Rydberg-type orbitals 
(6 a and 3 tt) are treated as high-level orbitals with M=1000, while some of the other 
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Rydberg-type orbitals (6 a and 10 vr) are treated as low-level orbitals with M=100. The 
results are listed in Table. It is shown that state-averaged ML-DMRG-SCF could give 
qualitatively correct and quantitatively reasonable descriptions of 2^A' and ?^A' states. 
The shift between these two states is determined to be 0.69 eV by ML-DMRG-SCF, 
close to that by the state-averaged DMRG-SCF calculation (0.75 eV). When the further 
DMRG calculation is implemented with the optimized orbitals from state-averaged ML- 
DMRG-SCF, we can get the results nearly the same as that by state-averaged DMRG- 
SCF. Additionally, we illustrate the calculated NOs of a” symmetry in Fig. [^(2^A', 3^A' 
states are belong to vr — vr* type excitation, thus only a" orbitals are involved). It can be 
seen that the NOs obtained from state-averaged ML-DMRG-SCF and those from state- 
averaged DMRG-SCF are similar with each other, whatever in the shapes or in the natural 
orbital occupation numbers (NOONs). 
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NOs by DMRG-SCF(10e,34o) 


^ 9 $ ## 

1.955 1.924 1.899 1.599 1.491 0.411 0.339 

-in 4|(H^ ^ ^ :44e i<S^ 

0.186 0.069 0.043 0.010 0.008 0.007 0.006 

9 ^ 

0.005 0.004 0.003 0.003 0.003 0.002 0.002 0.001 


NOs by ML-DMRG-SCF(10e,34o) 

1.968 1.944 1.929 1.583 1.539 0.412 0.384 





-?ii 4i!» itm *ii£ 


0.096 


0.050 


0.028 


0.006 


0.016 


0.008 


0.005 


^ 4# 4|||l|fr ^ 4|l» 

0.004 0.003 0.003 0.003 0.001 0.002 0.002 0.001 

Figure 7: The NOs as well as the NOONs of a" orbitals in indole molecule. Top: Or¬ 
bitals are optimized by state-averaged DMRG-SCF(10e,34o) with M=1000 (12a'-|-22a" 
orbitals). Bottom: Orbitals are optimized by state-averaged ML-DMRG-SCF(10e,34o) 
with M=1000 in high-level subspace (6a'-|-12a" orbitals), M=100 in low-level subspace 
(6a'-|-10a" orbitals). 
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Table 2: Excitation energies (in eV) of 2^A' and states (by state- 
averaged calculation for 5 states) in indole system and the energy shift 
values (in eV) between them. The geometry was optimized by BP86- 
D3/augTZ and the aug-cc-pVTZ basis set is used when implementing 
DMRG calculations. 



GASSCF“ 

DMRG- 

SCF'* 

ML-DMRG- 

SCF® 

DMRG / ML- 

DMRG-SCF'^ 

DMRG- 

SCF® 

exptl.-^ 

2M' 

4.67 

4.78 

4.87 

4.79 

4.76 

4.37 

3M' 

5.60 

5.58 

5.56 

5.53 

5.51 

4.77 

A 

0.93 

0.80 

0.69 

0.74 

0.75 

0.40 


CASSCF(10e,12o), Oa' and 12a" orbitals are chosen. 

DMRG-SCF(10e,18o) with M=1000, 6a' and 12a" orbitals are chosen. 
ML-DMRG-SCF(10e,34o), 6a' and 12a" orbitals are chosen as the high-level orbitals 
with M=1000, the other 6a' and 10a" orbitals are chosen as the low-level orbitals 
with M=100. The orbital rotations between different levels are not allowed. Typical 
time for per-SGF iteration is around 5h using HP ProLiant DL380p clusters (GPU 
: Xeon E5-2670x2, 16 cores). 

DMRG(10e,34o) calculation (M=1000) using optimized orbitals from c. 

® DMRG-SCF(10e,34o) calculation (M=1000). Typical time for per-SGF iteration is 
around llh using the same clusters as [c]. 


f ref. 115 
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4.4 Orbital ordering effect and entanglement analysis 


It is well-known that the orbital ordering would signihcantly affect the DMRG perfor¬ 
mances, and there were many informative suggestions 22,24 -26,34,35, ^ ^ 57,63,69 


on how to construct an optimal orbital ordering before implementing the computation¬ 
ally expensive DMRG calculations. When using canonical-type orbitals, one reasonable 
choice is to use the bonding/anti-bonding ordering, which has been recommended by 
many groups |^[25|[^ . However, in most cases it is not straightforward and it is also ex¬ 
tremely time-consuming to pick up the desired bonding/anti-bonding ordering manually, 
and as such, there are some other suggested automatic strategies, such as reverse Guthill- 


Mckee algorithm 22,24,116 , genetic algorithm 32 , Minimun Bandwidth by Perimeter 


Search 71 , and algorithms from graph theory 63,87 . Since the ML-DMRG strategy 
changes the orbital ordering when partitioning orbital spaces, it is of great importance to 
also investigate the ML-type orbital ordering effect on DMRG calculations. 


For illustrating the entanglement behavior and orbital ordering effect within ML- 
DMRG calculations, here we show our benchmark calculations of Gr 2 system at two 
geometrical configurations (inter-atomic distance 1.5A and 2.8 A). In our calculations cc- 
pVTZ basis is used with a 12e-in-42o active space. The NOs from GASPT2(12e,12o) are 
used and the partition of active spaces is as following: the valence orbitals (12 orbitals, 
i.e. 3d4s orbitals) are chosen as the high-level orbitals; 16 orbitals which is based on 
NOONs (mainly composed by 4p4d orbitals) are chosen as medium-level orbitals; other 
14 orbitals, which are also based on NOONs, are chosen as the low-level orbitals. Different 
Davidson convergence thresholds (10“^, 10“^, 10“®) are used for parallel comparisons. 

Fig.g and Fig. show the orbital ordering schemes and their partition of 1-body 
integrals as well as the entanglements (I’-orbital entropy and mutual information). Two 
facts can be easily found for the ML-DMRG orbital reordering strategy compared with the 
original ones: 1) the one-body integrals would no longer locate at the thick diagonal part 
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(Natural) 



Ag B3u B2u Big Blu B2g B3g Au 


Natural ordering 
in ML-DMRG spaces 
(ML-Natural) 



Low-level High-level Medium-level 


Figure 8: Schematic illustration of natural ordering and its ML expansion in Ct 2 (r=2.8 
A) systems with 12e-in-42o. Here one column represents one orbital and the height 
of the column reflects the NOON. Their one-body integrals (shown is squares) and the 
entanglements [I’-orbital entropy and mutual information (i.e. orbital interaction), shown 
in circles] are also illustrated. In the entanglements circles, the magnitude of red point 
reflects the single orbital entropy; the lines that link two orbitals represent the mutual 
information. The black solid line denotes the mutual information larger than 0.1, gray 
dashed line for 0.01, and green line for 0.001. The entanglements are calculated by 
Maquis (DMRG) using M=2000 states. 
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Figure 9: Same as Fig. schematic illustration of irrep blocks ordering, bonding/anti¬ 
bonding ordering in irrep block pairs, and their ML expansions in Ct 2 (r=2.8 A) system 
with 12e-in-42o. 
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but they will be split into several thin sub-bands; 2) the entanglements, both of the single¬ 
orbital entropy and mutual information, are no long well distributed in the whole orbital 
space, and they are concentrated in the high-level subspace. This may bring an advantage 
of quicker DMRG convergence in each subspace since the major entanglements are more 
localized than before. But it may also bring a disadvantage of slower DMRG convergence 


when accounting the unbalanced inter-subspace entanglements 26 . The comparison of 


DMRG convergence behaviors is illustrated in Fig. 10 for 2.8 A Gr-Gr distance and those 


for 1.5 A Gr-Gr distance in Fig. S2 110 . Interestingly, when loose or moderate Davidson 


threshold (10“'^ or 10“^) for wavefunction accuracy was used in DMRG calculations, the 
ML-type ordering may offer easier convergences to the lower energy than their original 
counterparts in the beginning, especially when the original orderings are well but not 
perfectly reordered (e.g. the irrep blocks ordering in our calculations). This can be 
ascribed to the above mentioned advantage of ML-type orderings in accounting intra-band 
interactions for moderately-accurate wavefunction calculations and also its shortcoming 
in describing the inter-band interactions. This provides a useful hint for future DMRG 
calculations: in the case of moderate accuracy requirement of the wavefunction, ML- 
type orbital ordering can be used as an efficient “warm-up” strategy to save substantial 
computational costs. 

When we further examine the efficiency difference between ML-DMRG calculations 
and DMRG calculations with hxed M values within the same hxed ML-bonding/anti¬ 


bonding orbital orderings, one can also hnd from the top panel of Fig. 11 that ML- 
DMRG calculations (with different ratio of preserved states shown in legend) do not 
have signihcant superiority over their counterparts. This can be understood from the 


discarded weight evolution pictures in the bottom panel of Fig. H, where the discarded 
weights in ML-DMRG calculations are generally much larger than DMRG calculations 
with hxed M values (blue line) in the case of the complicated Gr 2 system which is much 
different from the earlier N 2 system (see Fig. Sib 110|). This implies that ML-DMRG 
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calculations may be not efficient for very complicated correlation systems because such 
cases usually requires a large number of M at each orbital, whereas ML-type orderings 
can still be suggested for DMRG calculations with fixed M values for such complicated 
cases according to our above discussion. 





Kio’’ x10'’ 





Time (s) 


Figure 10: The convergence behaviors of DMRG calculations of Gr 2 with 2.8 A inter¬ 
atomic distance using different orbital orderings and different convergence thresholds. 
The Davidson threshold values are 10““^, 10“^, 10“® from left to the right. The M values 
were increased regularly: 50, 100, 200, 400, 800, 1200, 1600, and 2000. There are 4 
sweeps for each M value. The largest discarded weights during sweeps are also listed at 
the bottom. HP ProLiant DL380p clusters (GPU: Xeon E5-2670x2, 16 cores) are used in 
calculations. 
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Figure 11: The convergence behaviors of DMRG calculations of Cr 2 with 2.8 A inter¬ 
atomic distance using ratio in active subspaces and different convergence thresholds. The 


calculating condition is as the same as that shown in Fig. 10 with the bonding/anti¬ 
bonding ordering in irrep block pairs. Notice that only the energy region from -2086.7 to 
-2086.9 (in hartree) are shown for better presentations. 
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5 Conclusion 


In this paper, we build the DMRG algorithms in the ML-type framework and inves¬ 
tigate their efficiency and accuracy performances compared with their standard DMRG 
counterparts through computational examples of H 2 O, N 2 , Gr 2 and indole. Our major 
hndings can be outlined as follows: 

1. Ghemical intuition-based multi-level control of the large active space can be com¬ 
bined with DMRG calculations. A noticeable computational efficiency gain with a 
little price of energy inaccuracy (normally few mHartrees) can be obtained by ML- 
DMRG calculations in contrast to those standard DMRG calculations with hxed M 
values. 

2. The optimized orbitals obtained from ML-DMRG-SGF calculations are similar to 
those from the DMRG-GASSGF calculations. Gomparing to RASSGF, a larger 
high-level active subspace can be used in ML-DMRG-SGF. 

3. ML-type hierarchical orbital ordering could offer increased computational efficiency 
for DMRG calculations with medium accurate Davidson threshold, since the major 
intra-subspace entanglements are more localized within such a new ordering. How¬ 
ever, it will lose computational advantage for highly accurate Davidson threshold 
due to its inefficiency in describing the delocalized inter-subspace entanglements. 

The above findings concerning on the energetic accuracy and computational efficiency 
as well as the energy convergence behavior of ML-DMRG type calculations provide use¬ 
ful information for constructing more economic protocols for future large-scale DMRG 
calculations. For example, ML-DMRG-SGF could be taken as a pre-DMRG-GASSGF 
calculation in order to run a DMRG-GASSGF calculation more efficiently. The hierarchi¬ 
cal ML-bonding/anti-bonding orbital ordering may also be recommended for the DMRG 
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calculations in large and complicated active spaces. 
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